In [1]:

    

import os
import os.path as op
from numpy import arange
from numpy.random import permutation
import nibabel as nib
from surfer import Brain


    

In [7]:

    

%gui qt


    

In [4]:

    

subject_id = "fsaverage"
subjects_dir = '/Volumes/group/awagner/sgagnon/RM/data' # os.environ["SUBJECTS_DIR"]


    

In [29]:

    

coords = dict(green = [[-15,-66,54], #Hutchinson 2014 SPL (task>specific)
                      [-9,-72,57], #Hutchinson 2015 (decision uncertainty)
                      [-11,-57,52], #Fleck 2006 (low>high conf)
                      ],
                red = [[-30,-54,39], #Hutchinson 2014 (item only > cr)
                       [-39,-54,45], #Hutchinson 2015 (linear memory strength)
                       [-39,-55,36], #Wheeler/Buckner 2004 (K>N)
                       [-38,-76,28], #Daselaar 2006 (increase famil)
                       [-33,-60,36], #Yonelinas 2005 (increase famil)
                       [-3,-76,45], #Montaldi 2006 (increase famil)
                       [-39,-47,47], #Montaldi 2006 (increase famil)],
                       [-37,-59,42], #Sharot 2004 (K>N)
                       [-36,-80,33], #Sharot 2004 (K>N)
                       ],
                yellow = [[-39,-63,39], #Hutchinson 2014 (source>item)
                          [-50,-41,25], #Eldridge etal 2000 (R>K)
                          [-30,-79,43], #Fenker etal 2005 (R>K)
                          [-42,-72,39], #Henson etal 1999 (R>K)
                          [-43,-67,40], #Wheeler/Buckner 2004 (R>K)
                          [-53,-25,18], #Yonelinas 2005 (R>K)
                          [-42,-68,31], #Montaldi 2006 (R>K)
                          [-10,-68,29], #Sharot 2004 (R>K, R>N)
                          [-52,-59,24], #Sharot 2004 (R>K, R>N)
                         ],
                orange = [[-39,-30,45], #Vilberg/Rugg 2007 (famil (K>M, exclusive R>K))
                          [-30,-63,45]]) #Vilberg/Rugg 2007 (famil (K>M, exclusive R>K))])


    

In [30]:

    

hemi='lh'
surf='semi7'

b = Brain(subject_id, hemi, surf, background='white')

for color in coords.keys():
    for coord in coords[color]:
        b.add_foci(coord, map_surface="white", 
                   color=color, alpha=.8, scale_factor=.5)
    
b.show_view('parietal')


    

    
Out[30]:





((-119.99999999999999,
  59.999999999999993,
  419.26721191406256,
  array([ 0.,  0.,  0.])),
 49.10600000000001)

In [ ]: